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The emergence of drug-resistant bacteria is a severe problem in aquaculture. The ability of drug resistance genes 
to transfer from a bacterial cell to another is thought to be responsible for the wide dissemination of these genes in 
the aquaculture environment; however, little is known about the gene transfer mechanisms in marine bacteria. In this 
study, we show that a tetracycline-resistant strain of Photobacterium damselae subsp. damselae, isolated from seawater 
at a coastal aquaculture site in Japan, harbors a novel multiple drug resistance plasmid. This plasmid named pAQUl 
can be transferred to Escherichia coli by conjugation. Nucleotide sequencing showed that the plasmid was 204,052 
base pairs and contained 235 predicted coding sequences. Annotation showed that pAQUl did not have known repA, 
suggesting a new replicon, and contained seven drug resistance genes: WacARB-9-like, floR, mph(A)-\\ks, »?e/(A)-like, 
sul2, tet(M) and tet(B). The plasmid has a complete set of genes encoding the apparatus for the type IV secretion 
system with a unique duplication of traA. Phylogenetic analysis of the deduced amino acid sequence of relaxase 
encoded by tral in pAQUl demonstrated that the conjugative transfer system of the plasmid belongs to MOBhw, a 
sub-group of the MOB H plasmid family, closely related to the IncA/C type of plasmids and SXT/R391 widely distributed 
among species of Enterobacteriaceae and Vibrionaceae. Our data suggest that conjugative transfer is involved in 
horizontal gene transfer among marine bacteria and provide useful insights into the molecular basis for the dissemination 
of drug resistance genes among bacteria in the aquaculture environment. 
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MICROBES AND 
ENVIRONMENTS 



The emergence and spread of drug-resistant bacteria is a 
global public health problem. This problem is impacted by 
both human and non-human use of antimicrobial agents. For 
example, in aquaculture, antimicrobial agents are used 
prophylactically to protect cultured animals from infectious 
diseases and avoid economic losses. Use of antimicrobial 
agents leads to selection for drug-resistant bacteria (6). 

The occurrence of tetracycline-resistant bacteria in the 
sediment and seawater of an aquaculture site significantly 
increased after using oxytetracycline (38). Additionally, prior 
acquisition of drug resistance gene(s) by bacteria is a major 
factor for the generation and propagation of drug-resistant 
bacteria. Previously, the common tetracycline resistance 
gene, /e/(M), was found in several distantly related species 
of bacteria isolated from fish at an aquaculture site (23), as 
well as from the sediment and seawater at the same site (38); 
therefore, gene transfer among different species of bacteria 
appears to be involved in the acquisition of drug resistance 
gene(s) in the aquaculture setting. Previously, we found that 
the majority of the fe?(M)-positive isolates belonged to closely 
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related genera, Vibrio and Photobacterium, suggesting that 
bacteria from this group might serve as an important reservoir 
for tet(M) in the aquaculture environment (38), and tet(M) 
in some bacteria could be transferred to Escherichia coli and 
confer drug resistance (36). Among the bacteria from this 
group, P. damselae subsp. damselae is an indigenous marine 
bacterium (1, 46) known to be pathogenic in both fish and 
humans. In fish, it causes septicemia in a broad range of 
species, e.g., damselfish and rainbow trout (29, 40, 41). In 
humans, exposure to seawater or fish contaminated with P. 
damselae may cause fatal necrotizing fasciitis (2, 12, 17, 34, 
52, 59). Some clinical isolates of P. damselae subsp. damselae 
are drug-resistant, leading to obstacles for chemotherapy (59); 
therefore, drug resistance in aquatic bacteria may possibly 
pose a risk for human health either directly or indirectly. To 
evaluate this risk, determining the characteristics of drug 
resistance genes carried by P. damselae subsp. damselae and 
the transfer mechanism of these genes to human enteric 
bacteria is required. 

Here we determined the molecular basis for tetracycline 
resistance in a strain of Photobacterium damselae subsp. 
damselae isolated from seawater at a coastal aquaculture site 
in Japan. The results show that the strain has a transferable 
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plasmid named pAQUl, of approximately 200 kilo base 
pairs (kbp). Nucleotide sequencing showed that the plasmid 
contains seven drug resistance genes and a complete set of 
genes encoding the apparatus for the type IV secretion system 
(T4SS) with a unique duplication of traA. 

From the deduced amino acid sequence of the relaxase en- 
coded by tral, pAQUl was classified as a novel member of 
the MOB H plasmid family in sub-group MOB H i2- The phylo- 
genetic relationship of this plasmid compared to other known 
groups of plasmids and integrative conjugative elements 
widely distributed among the species of Enterobacteriaceae 
and Vibrionaceae is discussed. Many reports have shown 
the prevalence of drug resistance among fish pathogens (6, 
24, 33) and aquaculture environmental bacteria (15, 38), 
whereas there is little information concerning the mechanism 
of their emergence. Our data provide useful insights into 
the molecular basis for the dissemination of drug resistance 
genes among bacteria in the aquaculture environment and 
their clinical impact. 

Materials and Methods 

Bacterial strains 

P. damselae subsp. damselae strain 04Ya311 was isolated from 
seawater at a coastal aquaculture site in Japan and was used as the 
conjugation donor (38). Its identification was confirmed using a 
full-length 16S rRNA gene sequence and PCR detection of the 
ureC gene (39). E. coli K-12 strain W3 110 (whole genome sequence 
accession no. AP000091; 19) was obtained from the National 
BioResource Project (National Institute of Genetics, Japan) and used 
as the conjugation recipient. Rifampicin-resistant W3110Rif r , used 
as the recipient in the alternative transfer experiment, was artificially 
obtained by induction. The P. damselae subsp. damselae donor 
strain was cultured in brain heart infusion medium (BHI; BD, 
Franklin Lakes, NJ) with 2% NaCl and 20 ug mL 1 tetracycline at 
25°C. The transconjugant and recipient E. coli strains were cultured 
at 37°C in Luria Bertani medium (LB; BD) with and without 20 
ug mL 1 tetracycline, respectively. 

Filter mating 

The donor strain was cultured to the mid-log phase in BHI with 
2% (w/v) NaCl and 20 ug mL 1 tetracycline at 25°C. The recipient 
strain was cultured to the mid-log phase in LB at 37°C. Filter mating 
was performed on LB agar plates at 25°C without antibiotic selection 
for 12 h; and the transconjugants were selected as described 
previously (36). To see the transferability of pAQUl from these 
transconjugants to an alternative recipient, two transconjugants 
(TJW2 and TJW13) were used as donors and W31 10Rif r was used 
as the recipient. Filter mating was performed on LB agar plates at 
37°C for 6 h; colonies that grew on the LB plates containing 20 ug 
mL -1 tetracycline and 100 ug mL 1 rifampicin at 37°C were selected 
and transfer of tet(M) was checked by PCR (38). 

Antimicrobial susceptibility test 

Minimal inhibitory concentrations (MICs) of E. coli were 
determined using the broth dilution method according to the NCCLS 
guidelines (35). For P. damselae subsp. damselae, we used the same 
method with modifications of the salt concentration to 2% NaCl, 
incubation temperature (25°C) and culture period (48 h) appropriate 
for this microorganism (26). 

Pulsed field gel electrophoresis (PFGE) and hybridization 

The total DNA from the donor and recipient strains was 
fractionated using PFGE according to the SI nuclease method 
developed for plasmids larger than 100 kbp (3), and transferred onto 
a nylon membrane (GE Healthcare UK Ltd., Buckinghamshire, UK). 



The plasmid DNA was detected using Southern blot hybridization 
(50). A 637 base pairs (bp) PCR product of tet(M) was labeled 
with digoxigenin using a PCR DIG probe synthesis kit (Roche 
Diagnostics, Mannheim, Germany), and used as the hybridization 
probe. 

Sequencing and annotation 

Total DNA containing both of the genomic and plasmid DNAs 
was extracted from a selected transconjugant using the QuickGene 
DNA tissue kit (Fujifilm, Tokyo, Japan) according to the manufac- 
turer's protocol, and the purified DNA was sent to Takara Bio 
(Ohtsu, Japan) for sequencing. The whole sequence of the total 
DNA was determined by the pyro-sequencing method with a 
Genome Sequence FLX system using titanium chemistry (Roche 
Diagnostics). Approximately 5 ug of the total DNA was used for 
the standard Roche 454 shotgun sequencing method, and a set of 
data having 628,489 sequences with an average length of 367 
nucleotides was generated. The data were estimated to encompass 
25 times the entire DNA sequence. From the data set, the genomic 
sequence of E. coli W3110 was subtracted leaving the plasmid- 
derived sequences. The plasmid-derived sequences were then 
assembled de novo using the GS assembler (Newbler 2.0; Roche 
Diagnostics), and three sets of assembled contigs were generated. 
To fill the gaps between the contigs, PCR was performed using the 
total DNA template with LA Taq polymerase (Takara Bio), and 
with primers corresponding to the sequences close to the terminal 
ends of each contig. The nucleotide sequences of the PCR products 
were determined using the dideoxy chain termination method with 
a Big Dye terminator V 3. 1 cycle sequencing kit (Life Technologies 
Corporation, Carlsbad, CA). Fidelity of the final assembly of the 
circular plasmid was confirmed using the scanning PCR method 
with 26 sets of primers designed with web-based software GenoFrag 
(60). The primers were optimized for long-range PCR to amplify 
the 26 segments covering the entire length of the plasmid. 
Preliminary identification of protein coding regions was performed 
using the MetaGeneAnnotator (37). The data were compared using 
BLAST for sequences in the NCBI non-redundant protein database 
for length, identity, and coverage; the final annotation was manually 
performed using the In Silico Molecular Cloning Genomics Edition 
(In Silico Biology, Yokohama, Japan). 

Phylogenetic analysis 

A BLAST search was performed using the amino acid sequence 
of pAQUl -encoded relaxase (Tral) as the query sequence. The Tral 
sequence was aligned using CLUSTAL W and compared with those 
of 24 other known Tral in the MOB H family ( 1 6) (Table SI). Twenty- 
one highly conserved regions within the respective Tral sequences 
were joined using trimAL (8), generating 25 sets of approximately 
450 amino acid sequences (Fig. SI). Finally, we used the neighbor- 
joining (NJ) method with 1,000 bootstrap values, and similarities 
were compared using Molecular Evolutionary Genetics Analysis 5 
(56). 

Nucleotide sequence accession number 

The entire sequence of pAQUl is deposited in GenBank/DDBJ/ 
EMBL databases under accession no. AB571865. 

Results and Discussion 

Transfer and detection of a 200-kbp conjugative plasmid 
carrying tet(M) 

A previous study showed tet(M) was transferred from 
P. damselae subsp. damselae 04Ya311 to E. coli using 
conjugation with a transfer rate of (6.62±1 .61)x 10 3 (36); 
therefore, filter mating experiments were performed to obtain 
transconjugants. To examine the location of the tet(M) gene 
in the P. damselae subsp. damselae 04Ya311 genome, we 
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Fig. 1. Detection of plasmid pAQUl in the donor and a representa- 
tive transconjugant using PFGE (A) and Southern hybridization with 
the tet(M) probe (B). Lanes: M, DNA size standard (lambda ladder); 
1, representative transconjugant TJ311W2; 2, E. coli W3110; 3, P. 
damselae subsp. damselae 04Ya3 1 1 ; and 4, negative control strain of P. 
damselae subsp. damselae JCM8967. About 20 ng of DNA was loaded 
in each lane. 

performed PFGE and hybridization analysis using a tet(M) 
probe. The PFGE showed that an extra-chromosomal DNA 
of ~200-kbp was present in P. damselae subsp. damselae 
04Ya311, and the position of the tet(M) hybridization signal 



corresponded to this extra-chromosomal DNA (Fig. IB). 

PFGE and subsequent hybridization showed that extra- 
chromosomal DNA of -200-kbp was present in the 
tetracycline-resistant transconjugant E. coli; but not in the 
recipients, coli W3110 control (Fig. 1A, B). Twenty trans- 
conjugants were generated using independent filter mating 
experiments and were similarly examined, yielding the 
same data. Therefore, tet(M) appeared to be carried by a plas- 
mid whose size was ~200-kbp, unlike our previous report 
(36). In our previous study (36), we failed to detect this 
200-kbp plasmid because of the use of regular agarose gel 
electrophoresis, which cannot separate DNA fragments lager 
than approximately 20-kbp, instead of PFGE. 

General features of the pAQUl sequence 

To characterize the 200-kbp plasmid bearing tet(M), the 
nucleotide sequence of the total DNA extracted from the 
transconjugant E. coli was determined using the pyro- 
sequencing method. Because the entire sequence of E. coli 
W3110 was known, its sequence was subtracted from the 
whole sequence data. The three contigs after subtraction were 
shown to be continuous by using PCR, and further gap 
sequencing revealed the complete sequence of a circular 
plasmid, which was originally transferred from P. damselae 
subsp. damselae. The plasmid was named pAQUl; it 
contained 204,052 bp and 235 predicted coding sequences 
(CDSs). The gene organization and G+C contents of the 




Fig. 2. Circular map of pAQUl. CDSs outside and inside of the circle are coded clockwise and counter-clockwise, respectively. Putative 
functions of the products of the CDSs are indicated in color: red, replication, partition and termination; purple, DNA processing; blue, conjugative 
transfer; green, transposition or integration; orange, antibiotic resistance; yellow, other functions; and gray, unknown functions. The third circle 
indicates GC content where purple shows upper GC value above the center line and gray shows lower GC value below the center line. 
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Table 1. Classification of CDSs identified in pAQUl 



Putative function 


Similarity 
observed 


CDS no. and 
gene name 


Identity (%) 
(number of identical 
amino acids/total 
length) 


Query 
coverage 

(%) 


Protein similarity observed 


Replication, partition and 


Yes 


001 (repA) 


34 (76/226) 


62 


Putative replication protein 


LCI llHUd LIU 11 




037 (par A) 


76 (177/232) 


91 


Partition protein ParA 






038 (para) 


57 (233/407) 


100 


rp ■ , ■ 1 , ' -rr r~» 

1 ranscnptional repressor protein KorB 






191 (ter) 


60 (171/286) 


95 


DNA replication terminus site-binding protein 


DNA processing 


Yes 


033 
058 


96 (1054/1101) 

80 (579/725) 


99 
98 


Type I restriction enzyme 

r\\T * a * TTT 

DNA topoisomerase 111 






108 


76 (224/296) 


100 


DNA modification methylase 






113 


41 (55/134) 


97 


Type IIA topoisomerase 






128 


79 (158/201) 


72 


5 '-nucleotidase 






134 


60 (353/592) 


98 


DNA primase 






156 


A <-\ /01 / 1 n^> \ 

42 (81/193) 


94 


DNA topoisomerase I 






176 


C\f\ / O O / A 1 /\\ 

90 (388/430) 


100 


TT Tfc T A 1 ' , J T*\X 7" A 1 

RNA-directed DNA polymerase 






201 


71 (342/485) 


96 


DNA helicase 






229 


74 (306/412) 


99 


DNA methyltransterase 


Conjugative transfer 


Yes 


062 (tral) 


50 (528/1058) 


99 


Tral 






063 (traD) 


81 (504/621) 


99 


TraD 






065 (s043) 


60(116/194) 


100 


Hypothetical protein 






067 (traL) 


75 (157/209) 


100 


TraL 






068 (traE) 


68 (141/208) 


99 


TraE 






069 (traK) 


61 (192/314) 


99 


TraK 






070 (traE) 


57(251/443) 


98 


TraB 






071 (traV) 


63 (120/189) 


98 


TraV 






072 (traAl) 


70 (89/127) 


100 


TraA 






073 (traAI) 


67 (82/123) 


100 


TraA 






075 (dsbC) 


60 (122/203) 


85 


DsbC 






076 (traL) 


11 (629/814) 


99 


TraC 






078 (trnp) 


62 (106/170) 


99 


TrhF 






ATA / j fI7\ 

079 (traW) 


58 (242/417) 


98 


TraW 






081 (trau) 


82 (276/335) 


100 


TraU 






no 1 / j A 7\ 

082 (traN) 


66 (613/929) 


99 


TraN 






AA 1 / it Z" "7\ 

09 1 (s063) 


76 (156/205) 


91 


Hypothetical protein 






202 (trap) 


55 (l 74/3 1 8) 


88 


TraF 






i n"» / j t a 

203 (trati) 


70 (337/479) 


99 


TraH 






204 (traO) 


71 (860/1214) 


99 


TraG 


Transposition 


Yes 


050 


98 (365/373) 


100 


ISVhal (ISAsl family) 






056 


95 (229/240) 


100 


ISVba2 (IS5 family, IS903 group) 






157 


100 (312/312) 


100 


ISVsa5 (IS4 family, IS 10 group) 






158 


1 nn /"» i 1 /"> i 1 \ 

100 (311/311) 


89 


ion jC n st "> .c 'i \ 

lSSnlr9 (ln3 tamily) 






159 


AA /T>0/T>n\ 

99 (238/239) 


100 


irii; "> /Ton 1 .c "I \ 

ISVsa3 (1S91 tamily) 






163 


1 nn / A C\1 1 A C\1 \ 

100 (497/497) 


100 


1(117 "> /TOn 1 -C *1 \ 

ISVsa3 (1S91 tamily) 






1 7? 


qo H61/1 6M 


97 


ij v adj (lot lainiiy, laiu gioup^ 






173 


99 (236/237) 


100 


ISVsa5 (IS4 family, IS 10 group) 






174 


99 (586/587) 


97 


ISShfr9 (Tn3 family) 






178 


99 (401/402) 


100 


ISVsa5 (IS4 family, IS 10 group) 






187 


96 (298/312) 


100 


ISVsa5 (IS4 family, IS 10 group) 






189 


99 (233/234) 


93 


ISVsa5 (IS4 family, IS 10 group) 


Integration 


Yes 


024 


98 (309/315) 


100 


Integrase core domain 






025 


94 H 86/1 971 

y ^ 1 1 OU/ i J I } 


100 


iiiivgiujw lain ii y uiu lc in 






089 


78 (420/537) 


98 


Integrase family protein 






150 


60 (191/321) 


99 


Integrase family protein 


Antibiotic resistance 


Yes 


084 (fe/acARB-9-like) 


100 (223/223) 


100 


Beta lactamase class A (carbenicillinase) 






161 (floR) 


99 (402/404) 


100 


Florfenicol/chloramphenicol resistance protein 






165 (mpA(A)-like) 


63 (186/293) 


99 


Macrolide 2'-phosphotransferase Mph(A) 






166 (me/(A)-like) 


51 (208/407) 


99 


Macrolide efflux pump Mef(A) 






167 (sa/2) 


99 (270/271) 


100 


Dihydropteroate synthase 






177 (MM)) 


100 (639/639) 


100 


Ribosomal protection protein TetM 






181 (tet(B)) 


100 (401/401) 


100 


Tetracycline efflux protein TetB 
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Others Yes 002 


45 (640/1419) 


99 


Rhs family protein 


005 


52(154/296) 


97 


Rhs family protein 


019 


98 (235/240) 


92 


LabA-like protein 


029 


95 (483/508) 


100 


Type I restriction enzyme M protein 


031 


53 (225/428) 


99 


ATP-like protein 


034 


98 (252/258) 


100 


Metal-dependent hydrolase 


080 


40 (86/217) 


92 


Cyclic diguanylate phosphodiesterase domain- 
containing protein 


096 


78 (249/320) 


99 


Aerobic cobaltochelatase subunit Cobs 


102 


66 (400/605) 


100 


Von Willebrand factor type A domain-containing 
protein 


116 


58 (43/74) 


97 


Putative redox protein, regulator of disulfide 






bond formation 


131 


47 (131/278) 


92 


NgrC 


162 


100 (101/101) 


100 


LysR family transcriptional regulator 


170 


100 (237/237) 


100 


Mobilization protein B 


1 7? 

1 / J 




on 
oil 


Tn9/(5, transcriptional regulator (ORF9) 


1 01 

loz 


1UU (zUo/zUo ) 


1 fin 
1UU 


Tetracycline repressor protein TetR 


183 


100 (207/207) 


100 


ArsR family transcriptional regulator 


186 


99 (346/347) 


96 


Sodium/glutamate symport carrier protein 


188 


100 (58/58) 


100 


Sodium/glutamate symporter (fragment) 


206 


68 (48/71) 


86 


Ner-like DNA-binding protein 


209 


72(127/177) 


92 


Putative regulator protein 


222 


72 (228/316) 


98 


Periplasmic serine protease 


223 


57 (165/292) 


100 


DSBA-like thioredoxin domain protein 


Unknown Yes or No a 003, 004, 006-018, 


ND b 


ND 


ND 


020-023, 026-028, 030, 









032, 035, 036, 039-049, 
051-055, 059, 060, 061, 
064, 066, 074, 077, 083, 
085-088,090, 091-095, 
097-101, 103-107, 
109-112, 114, 115, 
117-127, 129, 130, 132, 
133, 135-149, 151-155, 
160, 164, 168, 169, 171, 
179, 180, 184, 185, 190, 
192-200, 205,207,208, 
210-221,224-228, 
230-235 

a Yes, similarity was observed with hypothetical proteins; No, similarity was not observed with known or hypothetical proteins. 
b ND, not determined. 



respective CDSs are summarized in Fig. 2. The CDSs 
were asymmetrically distributed: 182 CDSs were oriented 
clockwise, whereas only 53 CDSs were oriented counter- 
clockwise. When a BLAST search was performed using the 
deduced amino acid sequences of the CDSs, related proteins 
of known functions were found for 78 out of the 235 CDSs. 
Based on their predicted function, these CDSs were classified 
into eight groups (Table 1). For the other 157 CDSs, similarity 
to other proteins with known functions was not found (Table 
1). The average G+C content of pAQUl was 43.2%; however, 
there was a region of notably higher G+C contents at positions 
140,000 to 150,000 (Fig. 2) that contained transposons and 
the floR gene (see "Drug resistance gene" section). 

Genes potentially involved in replication of pAQUl 

Principally, plasmid-encoded Rep proteins have 3' — > 5' 
helicase activity and initiate plasmid replication through 
binding to specific sequences in the replication origin and 
facilitating the interaction of host proteins with the origin 
(55). The deduced amino acid sequence of CDS 001 showed 
33% (Table 1) identity to the putative replication protein 
(RepA) from Vibrio anguillarum and 24% identity to the 



Rep A protein of the well-characterized plasmid RSF1010 
from E. coli (31). Thus, although pAQUl appears to have a 
putative Rep protein, its sequence similarity with known 
RepA is low. Moreover, sequences similar to the primer sets 
conventionally used for PCR-based replicon typing (20) were 
not found to be effective for pAQUl. This suggests that 
pAQUl contains a novel replicon capable of replicating in 
both P. damselae subsp. damselae and E. coli. A previous 
report showed that unknown types of plasmids are distributed 
among marine environmental bacteria (11, 54). So far, little 
is known about the diversity of marine plasmids, although 
they seem to play an important role in the dissemination of 
resistance genes in the marine environment (15). 

Immediately downstream of CDS 001, there is an AT-rich 
region with five direct repeats and a set of inverted repeats 
(Fig. 3). The former is possibly an iteron, a region with AT- 
rich direct repeats that serves as the binding site for the 
plasmid-encoded Rep initiator protein (55). Such a replication 
origin containing an iteron is common to a large group of 
plasmids including the IncQ group plasmids (RSF1010, 
R1162, and R300B), PI, F, pSClOl, R6K, Rtsl, pColIV- 
K30, RK2, RP4, pCUl, pSa, and pPSlO, as well as in the 
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pAQU1_001_Rep 

ggcccatttg gtggtggcat ggcatcagag ctctccgaac ttcgtactgg taagtcctcc 

aagtaaj tcca a ttccttaaa ttaaqctaj: t gctttaacgc ct ttttcctt aaattaagct 

atatatacgg attaaggcgt ttcctcattt taagcttaga tataggtgtg tttcctaata 

ttaaqc ^ t atq ttatagagat qaac ttccta atattaaqc^t tcaa tttcct aatattaaqc 

tttagcttac agaaacaaag cttaaagtta ggaaat ttcc tqatattaag cta cctaaca 

cgctgcatga tcttatgatc ctatggatct atatatgatc agtagcgttc accgctgatg 

cggcttcgct gatcactcaa cttcgtttca ctgctatgag cagtaaaaca aagtcgagcg 

attcaattaa aacccgcttc catattagcc attgtttcat gtctgtatgc gcccgtatag 

ggcagttaag catgaagtta atctctttga ttaagttgat aatagagtcc gtgtagaaag 

taattaatgg actttagaaa ccagtggatc attttcatac taacttaatc gcaagagctt 

PAQU10Q2 

gt |ttgttagc acttatagtt cctagtgctg tttctgctca gtggcgtatc aatgccagtg" 

Fig. 3. DNA sequence of the potential replication origin region. 
CDSs are boxed. Arrows below the sequences indicate the 16 tol8-bp 
AT-rich direct repeats. Inverted repeats of 15 and 21 bases are indicated 
by dashed lines. 

lambda phage genome and in bacterial chromosomes (31, 
55). Among these, an IncP group plasmid, RK2, has the 
iteron-type replication origin downstream of its Rep gene 
(28, 43). The iteron-like structure of the AT-rich region in 
pAQUl and its position immediately downstream of the 
putative Rep protein gene suggests this region may corre- 
spond to the replication origin of pAQUl. 

The plasmid partitioning system consists of two proteins 
encoded by independent genes often known as par A and parB 
{korB), and a cis-acting centromere-like site, parS (14). The 
deduced amino acid sequences of CDSs 037 and 038 located 
30-kb downstream of CDS 001 were similar to ParA and 
ParB (KorB). The identities of the deduced amino acid 
sequences of CDSs 037 and 038 to those of ParA and 
ParB (KorB) proteins of the IncA/C plasmids were 75-76% 
and 54-58%, respectively, compared with those derived from 
P. damselae subsp. piscicida, (AB277723, AB277724), 
Yersinia ruckeri (CP000602), Yersinia pestis biovar 
Orientalis (CP000603), Salmonella enterica subsp. enterica 
serovar Newport (CP000604), E. coli (FJ621586, FJ621588, 
FJ705806), Salmonella enterica (AB277724), Aeromonas 
hydrophila (FJ705807), Xenorhabdus nematophila 
(FN667743), and A. salmonicida (CP000645). 

The amino acid sequence of the product of CDS 191 had 
60% identity to the DNA replication terminus site-binding 
protein (Tus) encoded by plasmids derived from Moritella 
sp. CDS 191 had 44^18% identity at the amino acid level to 
Tus encoded by the group of IncA/C plasmids mentioned 
above (except for A. salmonicida). A similar degree of identity 
was also found between the product of pAQUl CDS 191 and 
Tus encoded by an IncT plasmid derived from Proteus 
vulgaris (AP004237). 



Genes involved in conjugative transfer 

The deduced amino acid sequences of 18 CDSs (062, 063, 
067 through 073, 075, 076, 078, 079, 081, 082, and 202 
through 204) showed 50-82% identity to those of the 
apparatus of the type IV secretion system (Table 1). In 
addition, there were two CDSs (065 and 091) whose deduced 
amino acid sequences were similar to those of the proteins 
encoded by genes s043 and s063 found in SXT. SXT/R391 
is an integrative conjugative element (ICE) of less than 100 
kbp that has been detected in clinical and environmental 
isolates (4). s043 and s063 have been shown to be essential 
for conjugative transfer of SXT (58). pAQUl has all 20 CDSs 
encoding the proteins that appear to be necessary for the type 
IV secretion system-dependent transfer, suggesting that 
conjugative transfer of this plasmid is mediated by this 
system. Furthermore, observation of the transfer of pAQUl 
from transconjugants to the alternative recipient E. coli 
W3110Rif r (data was not shown) suggested that sufficient 
genes for this system were encoded on the plasmid. IncA/C 
plasmids also have both genes of the type IV secretion system 
proteins and genes similar to s043 and s063. Thus, IncA/C 
plasmids and SXT/R391 may have evolved from a common 
ancestral genetic element (58). Because pAQUl shares the 
organization of the 20 CDSs encoding proteins involved in 
conjugative transfer of the IncA/C plasmids and SXT/R391, 
pAQUl may have evolved from the same ancestor. 

Interestingly, pAQUl has a pair of traA genes located 
together, traAl and traA2, corresponding to CDSs 072 and 
073, respectively. The TraA protein, a TrbC/VirB2 homo- 
logue, is a subunit of bacterial pili serving as a pathway for 
DNA transfer from a donor cell to a recipient cell (27). The 
deduced amino acid sequences of TraAl and TraA2 proteins 
encoded by pAQUl are identical in length and 81% similar 
(101 residues/124). To incorporate TraA into mature pili, it 
is processed by removal of the 36-residue signal peptide and 
the 27-residue carboxyl terminal polypeptide (22). Amino 
acid sequences of the processed TraAl and TraA2 are highly 
conserved with only two amino acid differences (Fig. 4), 
suggesting they can both serve as functional subunits of pili. 
Phylogenetic analysis showed that pAQUl TraAl and TraA2 
belong to an independent cluster separated from other known 
TraA proteins (data not shown). It may be possible traAl 
and traA2 were generated by gene duplication within the 
plasmid. Another example of traA gene duplication is shown 
in the plasmid of Xenorhabdus nematophila ATCC 19061 
(accession no. FN667743). Although the biological signifi- 
cance of traA gene duplication is unknown, it is possible that 



Core 



TraA1_pAQU1 40 FDDIWTTITDWTQGTLGRIIAGAMILVGIVGGIARQSLTOFAMGIGGGVGLYNSPTWESIMTATLETADKVAPV 114 
TraA2_pAQU1 4 ° FDDIWTTITDWTQGTLGRIIAGEMILVGIVGGIARQSIMAFAMGIGGGVGLYHSPTWEAIMTATLESAEKIAPA 114 
TrbC_RP4 40 TGGSLPYESWLTNLRNSVTGPVAFALSIIGIWAGGVLIFGGELNAFFRTLIFLVLVMALLVGAQNVMSTFFGRG 114 




TraA1_pAQU1 mnvtstsnnkvmllfwvalalvisfmfanqahagtggte 

TraA2_pAQU1 MEITNKANNKANLLFHVGLVMIVGLMFANQAHAGTGGTE 

TrbC RP4 mttavpfrltmnrgilfylavffvlalalsahpamaseg 



AIQLSNGLGM 
AIQFSNGLGL 

AE I AALGNGALHQVQ VAAADAVRAVAAGRLA 



Fig. 4. 



Alignment of TraA homologues. Conserved amino acids are indicated by gray boxes. 
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this may lead to an increased level of TraA expression and 
pili formation. 

Plasmids encoding the proteins of the type IV secretion 
system have previously been found in pathogens for fish 
cultured in either seawater or freshwater, e.g. Aeromonas 
hydrophila, P. damselae subsp. piscicida, Y. ruckeri and A. 
salmonicida (13, 24, 45, 57). Our data show for the first time 
that bacteria carrying a conjugative plasmid could directly 
be isolated from seawater. Thus, conjugative transfer of the 
plasmid that carries multiple drug resistance genes may take 
place in the aquaculture environment and play an important 
role in the generation of multi-drug resistant pathogens. 

Drug resistance genes 

Deduced amino acid sequences of nine CDSs (084, 161, 
165, 166, 167, 177, 179, 180 and 181) were found to be 
homologous to those of known drug resistance genes. 

The nucleotide sequence of CDS 084 was identical to 
known but officially unnamed beta-lactamase genes 
(AB083415, AB453229) (Table 1) found in the plasmids of 
fish pathogen P. damselae subsp. piscicida strains isolated 
from cultured yellowtail (33). CDS 084 encodes a 224 amino- 
acid protein that has the serine-threonine-phenylalanine- 
lysine (STFK) tetrad active site characteristic of serine 
beta-lactamases (21). Further, its 55% similarity to a 
carbenicillin-hydrolyzing class A beta-lactamase (CARB-9) 
(44) and the presence of the arginine-serine-glycine (RSG) 
motif instead of the K-T/S-G motif of other class A beta- 
lactamases suggests that the product of CDS 084 is a 
novel carbenicillinase gene (WacARB-9-like) belonging to the 
CARB family (44). As described above, a previous work 
detected this gene in some strains of P. damselae subsp. 
piscicida, suggesting that these carbenicillinase genes may 
be widely distributed among Photobacterium in the marine 
environment. 

CDS 161 encodes a 405 amino-acid protein with 99% 
identity to FloR. FloR is a drug transporter of the major 
facilitator superfamily known as chloramphenicol and 
florfenicol resistance proteins. It is found in many Gram- 
negative bacteria including E. coli, K. pneumoniae, S. 
enterica, and V. cholerae. It is not uncommon for the 
transposase genes to be found upstream or downstream of 
the floR gene (51). DNA sequences upstream and downstream 
of CDS 161 are highly similar to the floR genes found in 
other plasmids and in the SXT element (51), appearing to 
contain putative transposase genes. CDS 161 is located in a 
region with a GC% much higher than the average of the 
entire plasmid, suggesting that the CDS may have been 
exogenously introduced to pAQUl by transposition. 



Based on the deduced amino acid sequence, CDSs 165 
and 166 showed 63% and 51% similarity to macrolide 2'- 
phosphotransferase Mph(A) and macrolide efflux pump 
Mef(A), respectively (9, 47), where the latter belongs to the 
major facilitator superfamily of efflux pumps. There is a 
standardizing rule for the nomenclature of macrolide resis- 
tance genes that the new gene must have <80% identity at 
the amino acid level with known macrolide resistance genes 
(47, 48); therefore, CDSs 165 and 166 are possibly new 
macrolide resistance genes, though their function needs to 
be confirmed to name them as new genes. CDS 167 encodes 
dihydropteroate synthase sul2 with 99% identity. This gene 
is found in many IncA/C plasmids (24, 57) and encodes 
dihydropteroate synthase with low affinity to sulfonamide, 
conferring resistance to the drug (53). The deduced amino 
acid sequences of CDSs 177 and 181 showed 100% similarity 
to the ribosomal protection protein and the efflux protein 
encoded by tet(M) and tet(B), respectively (48). The BLAST 
search identified the deduced amino acid sequences of CDSs 
179 and 180 as tetD and tetC (5), respectively, however, 
these were not characterized as tetracycline resistance 
determinants (18, 42), so we did not classify them as antibiotic 
resistance genes. Thus, many of the antibiotic resistance genes 
encoded on pAQUl have quite high similarity to other 
corresponding genes, while some appear to be new genes. 

Consistent with the annotation of the putative drug 
resistance genes, increases in MICs of the transconjugant 
were observed for ampicillin, carbenicillin, erythromycin, 
florfenicol and chloramphenicol, and plasmid-mediated 
transfer of drug resistance was confirmed (Table 2). Due to 
the inherent insusceptibility of the recipient strain of E. coli 
W3110, the transmission of resistance to sulfamethoxazole 
and sulfamonomethoxine could not be tested. Interestingly, 
all of the drug resistance genes found in pAQUl had a 
corresponding drug that has been permitted for use in 
aquaculture in Japan (32). Oxytetracycline and flumequine 
had indeed been administered in the aquaculture area where 
P. damselae subsp. damselae strain 04Ya3 1 1 was obtained 
(38). Most of the drug resistance genes are localized in 
between 140,000 and 160,000 positions and are flanked by 
transposase genes, apparently making an island region. The 
sequential transposition events and drug selection in the 
environment possibly played a role in constructing a multi- 
drug resistance island in the plasmid. Thus, plasmids 
possessed by bacteria in the aquaculture environment serve 
as an important reservoir for the various antibiotic resistance 
genes. The transfer of resistance genes among plasmids and/ 
or chromosomes of different bacteria may also occur in the 
environment. Impact of this ability to transfer among bacteria 



Table 2. Susceptibility of the donor strain of P. damselae subsp. damselae 04Ya311, the transconjugant TJ311W2, and E. coli W3110 to 
antimicrobial agents 



Strain name 


Species 






Minimal inhibitory concentration (ng/ml) a 






AMP 


CAR 


TET 


ERY 


CHL 


FLO 


SMXZ 


SMX 


04Ya311 


P. damselae subsp. damselae 


>128 


>128 


8 


128 


8 


16 


>512 


>512 


TJ311W2 


E. coli 


64 


>128 


16 


128 


8 


16 


>512 


>512 


W3110 


E. coli 


0.125 


1 


<0.125 


32 


0.5 


1 


>512 


>512 



a AMP, ampicillin; CAR, carbenicillin; TET, tetracycline; ERY, erythromycin; CHL, chloramphenicol; FLO, florfenicol; SMXZ, sulfamethoxazole; 
SMX, sulfamonomethoxine 
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TEH 



H3 



Name of plasmid or ICEs bacterial source (Inc group) 

R47S_Serratia marcescens (lncHI2) 
R27 _Salmonella fypft<{lncH1) 
Plasmidl Shewanella so ANA-3 



H11 



H12 



H2 



K;s Proteus vulgaris (IncT} 
pCA^ 1 Psntidcmanas resmovorans ;inr;P7) 
SXJ_Vibno cholerae 
R391 Procidentia rottgeri (IncJ) 
pAQU1 Phctobactenum damselae subsp damselae 
p\P V202_Yersinia pesVs biovar Onentalis sir P275 (IncA/C) 
pAsa4 Aeromonas salmomada subsp salmoniada A449 



clc genomic island Psoudomonus puddo 
PAG 2 Pseudomonas aeruginosa strain C 
l J AG -3_Pseucfomo/ias aeruginosa strain SG1 7M 
pKLC102_Pseuc/omonas aeruginosa 
PAG\-4_Pseudomonas aeruginosa strain C 
pMAQU02_Marv>7oftacfer aquaeolei VT8 
c hrom_ Neisseria gonorrhoeae 



Plasm id 1_Po/aromonas sp. JS666 
Plasm id ~\ Rhodoferax ferrireducens T1 1 8 
Plasmid2_Po/aromo/iassp. JS666 
pPNA01_Po/aromonas naphthalenivorans CJ2 

■ pBV\E02_Burkholderia vietnamiensis G4 
pHG1_Ra/sfo/iia eutropha H16 

Plasm \d2_Ralstonia metallidurans CH34 

■ pMOL28_Ra/sfon/a metaliidurans CH34 



CD 

£ 

CD 

o 



Fig. 5. Phylogenetic tree of relaxases (Tral) from plasmids and SXT/R391 that belong to the MOBn family. 





3H motif 




HD motif 




| pAQU1 1 


PVIGRYAAFVHLLPASEAHHHRGAGGLFRHGLE 


143 


LLHDVGKPLSD 


192 


plP1202 


PVIQRYAAFVHLLPASESHHHRGAGGLFRHGLE 


145 


LLHDVGKPLSD 


194 


pAsa4 


PVIHRYAAFVHLLPASEAHHHRGAGGLFRHGLE 


138 


LLHDVGKPLAD 


187 


SXT 


PCIDHLAAYVHLLPASEHHHHSGAGGLLRHSLE 


133 


LFHDIGKPVSD 


182 


R27 


KPLIEYIRWIHLLPASENHHHNGIGGLLSHSLE 


178 


LVHDAGKVYDL 


227 


R478 


KPLVEYVRWIHLLPASENHHHNGIGGLLSHSLE 


180 


LVHDAGKVYDL 


229 


pCAR1 


PILVNFAELAHMLPASENHHHSGPGGLLRHSLE 


212 


LLHDAGKPLTD 


261 


Rts 


PMLTNFANFVHLLPASEFHHHRAQGGLLRHTLE 


149 


VMHDIGKPISD 


198 


Plasmid1_Shesp 


DWANFLSYAHLLPASENHHHAAVGGLARHSLE 


171 


LVHDIGKALHD 


220 




APLERYAELVQRFPASEAHHHAYPGGMLDHGLE 


123 


LLHDVGKLAVD 


172 


PAGI-2°C 


APLERYAELVQAFPASESHHHAYPGGMLDHGLE 


123 


LLHDIGKIAVD 


172 


PAGI-3 SG 


APLERYASLVQQFPASESHHHAYSGGMLDHGLE 


123 


LLHDIGKVAVD 


172 


PAGI-4 C 


RPLANYAELVQQLPASENHHHAHPGGMIDHGLE 


89 


LAHDIGKI WD 


138 


pKLC102 


RPLANYAELVQQLPASENHHHAHPGGMIDHGLE 


120 


LAHDLGKIWD 


169 


pMAQU02 


YAIEQFAELVQEVPASEIHHHSESGGLIDHTLE 


127 


LSHDLGKIVTD 


176 


Chrom_Neigo 


PFLEKYIAFVQRLPASESHHHAGDGGLVRHTLD 


155 


LLHDVGKTVTG 


204 



Fig. 6. Conserved sequence motifs in relaxases (Tral) grouped in the MOBn family. Both of the 3H and HD hydrolase motifs indicated by gray 
boxes are conserved in the relaxase encoded by pAQUl. 



needs to be investigated in relation to the emergence of 
resistance in fish pathogens (23) or aquatic environmental 
bacteria (15, 38). 

Phylogenetic position of pAQUl 

To characterize and classify the plasmids derived from 
human clinical isolates, replicon typing using DNA hybrid- 
ization or PCR can be used (10). However, replicon typing 
tends to be difficult to apply to plasmids isolated from marine 
bacteria because their classification does not fit the plasmids 
derived from human clinical isolates (11, 54). Recently, 
comparison of the primary structures of tral gene products 
(relaxase) was successfully used for phylogenetic analysis 
and classification of the conjugative transfer system in six 
MOB families. Therefore, comparison of relaxases may be 
useful for plasmid classification (16). Relaxase is an essential 
enzyme for the type IV secretion system and cleaves at a 
specific site called oriT when it is transferred to another 
bacterial cell by conjugation. Because the replicon type for 
pAQUl could not be determined as mentioned above, 
phylogenetic analysis of pAQUl was performed based on 
the deduced amino acid sequence of relaxase it encodes. 

The data indicate that pAQU 1 is a member of the MOBh 
family and grouped in MOBhi2 together with IncJ, IncT, 
IncP7, and IncA/C (Fig. 5). The relaxase encoded by pAQUl 
has both 3H and HD motifs [(H/Q)-X 2 -PASE-X-HHH-X 3 - 



GG-X 3 -H-X-L and (L/V)-X-HD-(A/V/L/I)-GK, respectively] 
known to be shared by relaxases of the MOBhi and MOBm 
clades (Fig. 6). Among the MOBhw plasmids, pAQUl is 
most closely related to IncA/C plasmid and SXT/R391. In 
addition to the organizations of the 20 CDSs encoding the 
type IV secretion system proteins, the primary structure of 
the relaxase also demonstrated a close evolutionary relation- 
ship of pAQUl with the Inc A/C plasmid and SXT/R391; 
however, the DNA sequence of the tral gene of pAQUl is 
clearly different from those of the IncA/C plasmids and SXT/ 
R39 1 . Recently, evidence suggested that the IncA/C plasmids 
(7, 25, 30, 57) and SXT/R391 (4) found in various ecological 
niches may be responsible for the generation of multi-drug 
resistant bacteria. Similarly, pAQUl, also classified in 
MOBhi2, may play a similar role in the aquaculture 
environment, as well as in clinical settings. To examine this 
possibility, further investigations concerning the distribution 
of the 'pAQUl -type' plasmid will be necessary. 
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